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ABSTRACT 


There are many factors, such as aircraft configuration 
and weight, winds aloft, airspeeds flown, altitude, distance, 
etc., which affect fuel consumption in turbojet aircraft. 

For any given combination of these factors a flight path can 
be determined that will result in the least fuel consumed for 
a ground distance covered, Under divert conditions from air= 
craft carriers at sea to fields ashore the choice of the 
optimal flight patn is critical, The many possible combi- 
nations of factors lead to the adoption of computer flight 
planning. Pilots can avail themselves of computer solutions 
during flight planning and briefing sessions, and after 
take-off can receive further information via UHF radio. 
Typical flight handbooks display fuel flow data, etc. in 

Such a manner that the pilot must "guesstimate" entry 
parameters such as average horizontal weight, or weight prior 
to descent. Several iterative prccedures are developed that 
provide exact solutions to these important figures. MThus 

Che computer flight planning system will provide more 
accurate solutions, and free the pilot from this chore so 


that he may better spend his time briefing tactics. 
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CHAPTER I 
PURPOSE aOr Th SslLuDy 


Naval Aviation has been described as hours and hours 
of boredom interrupted by moments of sheer terror. The 
Simple word " BINGO " provides for some of these anxious 
moments, Bingo is well known to everyone as a friendly 
game of chance, but the word carries a special meaning to a 
U.S. Navy carrier pilot. It causes concern both ashore and 
at sea, and especially in the cockpit of at least one air- 
plane that is flying circles around some aircraft carrier 
somewhere at sea, 

"Fortress 501, this is Atlas tower, Bingo 100°/125 
miles, Change to departure control, 316.6, for radar 
mee vor, " 

"Atlas tower, this is Fortress 501. Roger, switching 
pio. o. " 

This conversation might have occurred west of San 
Clemente Island off the coast of Southern California. To 
Che pilot of Fortress 501 the message is unmistakable. de 
is to turn to 100° magnetic and land at the field 125 miles 
away, rather than attempt to land aboard Atlas, U.S.S. 
INTREPID, CVA-ll. 

In most cases an aircraft will receive a BINGO because 
of a flight deck accident, or because the pilot is having 


difficulty landing aboard the carrier due to bad weather 


and/or a badly pitching flight deck. The high fuel con- 
Sumption of jet aircraft at low altitudes precludes having 
the pilot hold until the deck is clear, or in the latter 
case attempting more landings aboard, 

When the pilot of Fortress 501 switches to 316.6, 
he will receive a radar vector to the nearest suitable 
landing field, the weather at the fieid, and the approach 
control and tower frequencies, But he will have to detere- 
mine his own flight profile; tnat is, how high to climb, 
What airspeeds to fly, and at what distance from his 
destination to commence an idle descent. Since the air- 
Craft is usually low on fuel when the pilot commences the 
BINGO, the pilot attempts to fly the flight profile that 
minimizes fuel burned, There are no other constraints gon 
the pilot, for he is only interested in running out of 
ocean before he runs out of fuel. 

The factors that affect fuel consumption in a jet 
aircraft are altitude, aircraft weight, and of course, the 
air distance travelled, Jet aircraft burn enormous amounts 
of fuel at low altitudes and during the climb to higher 
altitudes, During an idle descent from altitude the fuel 
flow is only one-sixth that of the climb portion. A 
heavier aircraft burns more than a lighter one, and a head- 
wind increases air distance and thus increases fuel burned. 
With these factors in mind, the pilot must decide what com- 
bination of climb, cruise, and descent flight paths will 


Cake nim to his destination with the Least fuel consumed. 


A wrong choice could result in the loss of an aircraft and 
perhaps a pilot. A lesser error might result in a safe 
landing, but several hundred pounds of jet fuel will have 
been wasted, This latter error might also result in some 
premature gray hairs on an intrepid naval aviator. 

Since a determination had to be made at what fuel 
weight to BINGO the pilot, the knowledge that all pilots 
would fly an optimal (minimum fuel) flight path would allow 
the BINGO to be delayed for perhaps one more pass at the 
flight deck, Also, any “gravy” could be reduced from the 
bingo weight. This lower weight would increase operational 
readiness by keeping all aircraft aboard the carrier where 
they belong, and where all commanders want them. If an 
aircraft is sent to the beach, the carrier must remain in 
the area to send messages concerning ship's position and 
overhead times so that the stray bird can come home to 
roost. Since carrier skippers like to hide their ships 
out at sea, there are obvious tactical advantages in 
keeping all aircraft aboard. 

An ideal way to ensure that all pilots fly an optimal 
flight path for all distances, aircraft weight, and wind 
conditions would be by computer flight planning. Thus when 
Fortress 501 contacted departure control on 316.6 he might 
hear, 

"Roger, Fortress 501, Climb at 285 knots to 19,000 


feet. Cruise at 260 knots indicated, Commence idle descent 


when 44 miles from destination fix." Departure control could 
then follow with the standard weather and radio aids infor- 
mation. 

This thesis considers only the short BINGO problem in 
its development of a computer flight planning program. This 
applies to a small, but important segment of aviation. The 
obvious extension is to use this type of flight planning for 
all long distance cross~country flights. The only modifi- 
cation necessary is to incorporate several horizontal legs 
at altitude sandwiched between a climb and a descent leg, 
rather than just one as does the BINGO program. The cross-~ 
country program would also require a change of optimal 
flight altitude to a higher one when the aircraft weight 
decreases enough to Warrant a climb in order to remain 


optimal. 
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CHAPTER II 
REPO OF Tk STUDY 


Tne BINGO problem may be stated mathematically as a 


Classic minimization problem: 


min (Wy 7 W,) 


Se Us h <40 x 10, in even thousands 


xX. of x < A, - 4%, 


Ne > B where B= aircraft empty weignt 
Wy No 
| | 
| | 
fn 
| | 
| | 
Wo , | , 4 
Ae x A 
Figure 1 


Definitions and abbreviations: 


Wo = aircraft weight at start of flight 


W = weight at end of climb portion 
No = weight at end of horizontal portion 
q. = weight at destination 
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eon 


rT MN 
= 


CAS 


TAS 


Ph 


GSPD 


7 f 


D1 


altitude in thousands of feet 

ground covered in climb portion 
ground covered in horizontal portion 
ground covered during idle descent 
total distance to destination 

time 


Calibrated airspeed, velocity of the air over 
the pitot tube 


true airspeed, CAS corrected to sea level 


iy 
TAS = — » where d= aaa 
d Po 


density of the air at altitude h 
density of air at sea level 


ground speed, aircraft speed relative to the 
ground (TAS corrected for wind) 


fuel burned in pounds 
rate of fuel ,fiow in: pounds per nautical mite 


drag index for the varticular aircraft con- 
figuration 


A flight profile depicted by rigure (1) is used because 


of standard operating procedure, and also because the air- 
craft instrumentation is such that it is the only profile 
that a pilot can fly accurately. An arcing path such as a 


semi-circle or a cycloid is impossible. 


In order to fly a flight profile in the form of a 


cycloid or a semi-circle a pilot would have to rely on his 
VSI (vertical speed indicator) to accurately control his 


rate of climb throughout the entire flight, This instrument 


de 


is inaccurate at best since it is a pressure instrument and 
fluctuates with slight changes in pressure. ‘The VSI also 
tends to lag the true pressure changes. A climb, level, 
descent flight profile can be accurately flown because the 
pilot need only control his airspeed and altitude, Both 
these cockpit indications are excellent in all aircraft, and 
a pilot gets acquainted with them from his first step into a 
cockpit. 

The objective function, Wo) = el is the fuel consumed 
in travelling the distance X,, Since the only weight loss 
Will be Glee fuel burned; The constraint W, > B is obvious, 
If the destination cannot be reached, even by flying an 
Optimal flight profile, then the aircraft should rendezvous 
With a tanker aircraft if one is available, or eise continue 
landing attempts aboard the carrier. It is not obvious that 
XK, + A < K ~ X, 1s a necessary constraint. A possible 
flight path could be to continue horizontally until the 
destination is reached, and then commence a circling descent 
such that X 


d 
descent at the proper airspeed for Xqg will use less fuel 


= 0. Navy flight tests show that an idle 


than any other possible flight path covering the distance 


Xx This type of idle descent is standard operating proce- 


a 
dure for jet aircraft. 

The altitude, h, must necessarily be less than or 
equal to the service ceiling of the aircraft being flown. 


For this study h < 40,000 ft., the service ceiling for the 


A4C, is used. Even thousands of feet are used, resulting 


in 40 possible flight profiles. Using every 500 feet will 
result in 80 possible profiles, and using every 250 feet 
Will result in 160, This study uses even thousands merely 
for ease of presentation. Only a minor change is required 
to solve for any desired number of possible profiles. 

The following assumptions were made concerning the 
filieht path of Figure I; 

1, the aircraft makes an instantaneous transition 
from the climb to the level flight attitude, and 
from the level to the descent, 

2, all airspeed changes are instantaneous, 

3. the aircraft can remain on track despite any 
cross-wind. 

4%, wind information is known, and there are no up- 
drafts or downdrafts. 

Performance data such as fuel flow and airspeed were 
obtained from NATOPS Flight Manual? for the A4YC aircraft. 
Although only the A4C is considered here, similar data can 
be used for other aircraft in the inventory. All perform=- 
ance data was gathered by U.S. Navy flight tests. The 
peculiar form of the graphs (Figure 6-13) makes it diffi- 
cult to form mathematical functions for such values as 
Maximum range airspeed, This precludes formulating the 
problem as a standard Lagrange minimization probien. 

The computer program developed is one that calculates 
the fuel burned for all altitudes up to 40,000 feet, the 


service ceiling of the A4YC, and chooses that flight profile 


that results in the least fuel burned. The first altitude 
tried is level at 1000 feet for the entire distance, Ko 
The remaining thirty-nine are flown as follows. 
(1) a 100% power climb to h-thousand feet 
(2) level at n-thousand feet at the max-range air- 
speed, (that airspeed that gives the most miles 
per pound of fuel) 
(3) an idle descent at that airspeed that covers 
Xs With tne least fuel burned 
Aircraft performance data is presented as follows. 
CLIMB PORTION: 


Iss 


I 


F (Wos ede, 
CAS =C (Ios Dao) 


time 


fi 


t (los h, DI) 


HORIZONTAL PORTION: 
ay +a 





f =f ( ne » Ny DE 
2 
ue Xn = 4 we | 1 
Pp = 5 Wy o- a, (1) 
Ny kW 
F = F(X; wig eae) 2) 
W + iy 
CAS = O(H-& » hy DI) 


DESCENT PORTION: 


CAS = C(wW,, DI) (3) 
X; = X(W5, h, DI) (4) 
time = t(h, DI) (5) 
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Future functional notation will omit the drag index, DI. 
It will be held constant at DI = 50, the value for an air- 
craft with no external stores. This is almost always the 
case in a BINGO situation, 

The first flight profile tried is level at h = 1000 
feet. Since this differs from the others only insofar as 
there are no climb or descent portions, a typical climb, 
level, descent profile will adequately explain both possi- 
bilities, For the climb to any alvitvude, Wo and h are known, 
From Figures 6, 7, and 8, F, CAS, and the time to climb 
are easily determined, CAS is converted to TAS, and tne 
predicted wind information (averaged over h=-thousand feet) 
converts TAS to GSPD. Wo - = Was and GSPD x time = Xo 
Bach 5000 feet the CAS is changed to the optimal climb 
speed for the next 5000 foot portion, If h > 5 the climb 
portion of the flight results in ([h/5]+ 1) iterations, 


At the start of the horizontal portion, W., h, and 


1 
(X, ~ x are known. Figures 9 and 10 however, require the 
average weight during the horizontal leg, (Wh + Wi )/es as 

an entry parameter, Recall that the fuel burned during 

the horizontal leg is an implicit function of No and Wa 
From equation 2, 


| ype ame hi eka Ae ote 


ol ve 
The procedure for using the graph in Figure 10 is to 


enter with the average weight, proceed horizontally to the 


16 


altitude h, read vertically downward to the DI = 50 line, 
and horizontally to the left to read f, the fuel flow, 
Three different methods were tried to solve the im- 


DIACa LelTUReChIOn FOreyeNe average Weilznit, 


METHOD I 
Approximate the altitude lines in Figure 10 with 


straight lines as shown in Figure <c. 


DI = 50 





Pieuwre 2 


1? 


From geometry we obtain: 


es f 
tan @ = ———=—_—— = ——_—_——— 
CW, — dO) ,000) Ny i Ie Wo 
eae 
a CW = +) J. Be 6 


Knowing t, we solve for fuel flow, f, using 
f =f cl ee Ea 
O 


Va of W 
a 5) tan @ tana 


o 


- 
| 

an 
4 
= 


(0) 





a ae ss 
re) L i a 


where N. is defined in Figure 2, 


Solving this quadratic for W. will allow a good approximation 


N+ | 
to the average weight, (422). This approximation is 
believed to be quite accurate, since the altitude lines are 


fairly straight. An exact approach is given in metnod two, 


HETAOD ie: 

Allowing for the fact that the altitude lines are non- 
linear, a more precise but lengthy method for solving for 
Che average weight is the method of successive approximations, 
Referring again to tne graph of Figure 10, we see that the 


average Weight is a number such that, if it is used to enter 


ViSmeaoke si beiroduces, a Fuel used. Fs such that W, ~ P= W,. 


W 


Wa bow 
nen Ww, is combined with W, to get the average (=~), this 


Zz 1 
average Will exactly equal the average we used to enter the 


table. Because of the construction of the table, the true 
average is the only entry number that will give the average 
back again. 

The first approximation to the average is We The 
aircraft continuously loses weight, thus we know that this 


approximation is too high. Entering the table with sr yields 


an FF, This yields an average 1? # wo Decrease 


A vo 


a, by some Aj >O. Entering the table with W,7 4 gives a 


Form another average: 


W + (W.=- F 


) 
a lag 


Continue in this manner until 


Ny a CW _ Sate 
a = Ue 
2 1 
This method is easily suited for computations on a digital 
computer, but it proved too lengthy. The greater the accur- 
acy desired, the smaller A must be, But a small A requires 


a large number of iterations, This metnod Was abandoned 


Hom tie third and final method, 


MaTHOD Lis 
Wy + We | 
Recall from Figure i0, dl aaa wij. The graph 


shows that the fuel burned is directly proportional to the 


aircraft weight As a first approximation to the true 


average weight, use Ne This-“Vs'-ctoo-lareée,. tnerel ores 


results in a fuel burned, Ps that is too large. 


f= F( AV h) 


1? 
From equation 6 
pecan tee 2 ase 
1 fo+ (Wy - AVS) tan @ tan a 
Let 
x = D a constant 


a" of W tan @ tana = B a constant 


tan @ tana = C a constant 


The first approximation to the fuel used on the horizontal 


leg is; 


D 
ee where X. = AV 
fy B+ Cx, 7 T =a (7) 


Selig lexs. Pr is too large, as a next approximation to the average 


weight decrease AV, by coe Thus in equation 7, multiply by 


- 3, and add dW, to both sides, we get 


=D 
2 1 2(B+ CX, ) ca 


= 
< 
ii 
_ 
nok! 
) 


A 
pee secs Cope t7 
20 "2 Becx, 1 
(where A = = =) 


Biber ine Preume 10) wi tn 4, yields 


_ D 
& B+ CX, 
and 
AV = v8 = a eg W 
3 ieee ere i 


2 


This iterative process yields a sequence {An} . 


X = W 


1 L 
x A 
2 B+CX, 1 
A 

A = t wW 
2 B+ CX, 1 
i ee 
nN ; % 

B+CX, , i 

wnere 
eS ealsie 
B= .34 
ay, 

C= 6. x 10 
Wose 1S Xe 10° 


aw 


Solving for tne first few points yields a nest of closed 


intervals as shown in Figure 3. 


Pie uses 


At this point it can be seen why (A, ~ BF, ) was used 
i 


as the second approximation to the average weight. Since 





we know that x is well above the average, it is desired toa 
decrease it by a quantity large enough to drive thewseead 
approximation below the average. A number like .8 or .9 
could have been used to multiply Fs but the use of $ is 
sufficient and causes the sequence to converge iaggea than | 
numbers like .8 or .9. 

The manner in which the successive points alveriauc 
led to applying the theory of continued fractions in an 
avvemp’ to show that the sequence ek converges. Hall and 
Knient’*- snow that each successive convergent of a con- 
tinued fraction is alternately less than and greater than 


the true value of the continued fraction, 


Rewriting equation 8 yields 


ee 


Back~substituting for x, ] by X93 for X, 2 by x. 3) 


etc. yields; 


Te ae A 
a ae a 
B+cC (NV, + ——————-——___ 
il Q 


B+C (Ww, + ————— 
i oe eee 
oe 
gee he. Sees 
B+ CX 
O 


mnesform of thevelassic Continued fraction 1s obvious. 
Factoring out the A in tne numerator, and multiplying in 


each convergent by C yields: 


X™ = Lim Ky = NJ +A a 
aaa C20)... ——— 
(Be OWs jy eee 
= (B+ CW, ) + 


Dividing numerator and denominator of each successive 
meacbion by (CA) 


ce 1 


To se 
bee ae al 


ae er ne ec ee - core pemmecene 


CA (B-- Ci ) + 


dl B-F Cul, 
eee Ee 
Ce 0 


In the continued fraction in the brackets, let 


a= B-+ CW >0 
a 
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Tnen 








5* 7 1 
oe i? 
D+ + 
a 
lens 
As we let n-eo, replace the ard ch | eve. 
o* t 
Be 
b + @ 
ena ae 
a(b 08) + 1 a(b+O*) +1 
aren 
b + 6 
Solving the resulting quadratic yields: 
5 oF 
5* . -abs (a@b* + bab)? 
2a, 
Thus we have that 


The minus sign is selected because ee Gy 


This third method gives the average weight for any 
given value of W, x? Sanda. 


converges in 6 to 8 iterations. whereas the second method 


< a 
-ab ~ (a@p* + 4ab)® 
1 2a, 


Trequitgn Pere gaan 50. 


Z'4 


i: 


For most values, the sequen@ 


(9) 


convergents by 











a | 
| 


1 


Al 
| q 











It is interesting to note that in the continued frac- 
tion, equation 9, for real but unequal odd and even con- 
vergents, Van Vleck's tneorem can be applied. For aor b 


imaginary, Stieltjes theorem applies. Proofs of these 


[3 


theorems are given in dall o Having solved the implicit 
W of N 
mounmetion for ———s Figure 9 yields CAS. It is converted 


to TAS, and the wind at altitude h changes TAS to GAPD 


time = xy / GSPD 
GAtdastanece. = 1S = time x TAS 


Ny of No 
Figure 10 gives F as a function of x, and =o 


Since optimal CAS is a function of the aircraft weight, 
as the weight decreases during the horizontal leg the CAS 
should be changed to remain optimal. In this program, when 
the aircraft weight decreases by 500 pounds, CAS is recon- 
puted. This results in breaking up x, into segments of 
length X = 500 xf. This results in ([X,/X, J + 1) iter- 
ations of the horizontal leg computations. - 

AG this point it should be noted that in order to 
solve for the fuel burned in the horizontal leg by either 


of the three methods, we require X At the start of the 


as 
leg we know only X, and x. Since X = xX, _ (X,+X,) the 
Veluwe Or xX. must be determined prior to the start of 


qd 
the horizontal leg. From equation 4 recall that 


Ka 7 CW 3 h). The problem resolves to this: in order 


to solve for We, we need X But we need W. before we can 


d° ed 
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soive for Kae Ancther method of successive approximations 
is developed that results in a recurrent sequence that is 
shown to converge. Once it is shown that we can find x5 
(and consequently XK) at the start of the horizontal leg, 
then all previous procedures are justified, 

The air distance covered during tne descent is 
og = X (h, DI, wind). From equation 3, CAS = €¢ (Wes Dit 4 
HOM -cauau vom oe b= =the, Oh) 


Prom Figure ll 


245 a } 


CAS 
190 3 | 
| 
| 
12,000 20,0008 
W 
Pa 
Figure 4 
CAS = 190 + tan gp (WwW 12000 ) 


Let (190 = 12000 tan gj =a 


{ 





and 


LAS = i meCoor ke “eo °CAS Where k = i. 
d O a 


Aad Ee a 


c ae 2 = ae TT a a 


d “Tas  GSPD 


if 


Cime in descent 


=U 
x a XeGorD 


Geeo = 204s 23nd et Kelas 
where x, is a constant to correct TAS for wind, 


GSPD 


ii 


ky (x CAS ) 


k (x La + bu, J) 


G3PD 1 


i 


Since 


X. = t_ X GSPD 
4 Je 


d 

hte ak _k + bk_k W 

a d ( ime) iG 2? 
Let 

ou = 15) 

x — 1. oe 3. 10 

A Bi, (10) 
Begin the iterative process by letting Wd = ane Ries 
large value will yield an Ko toat 18 LOO) large. ~ since 
x = x = (X an Ka )s the resulting x Will be too small. 


Fuel burned on the horizontal leg is directly proportional 


to the distance, therefore f Will be smaller than the true 


Zi 


= J - : his value is st t | 
Pe 1, B This value is still too Saree 


is closer to the true value of US Chan was the first approx 


value. ( 


mat Lon % Wy ° 


From equation 6 we have 


| Ay i 
qj es | 
é f +wW, tandtana = fan @ tama «/ Tg | | 
fe) i me | 
Let 
f + W. tan @ tan@®@ = B | 
QO a 
and 
- tan 0tan @ _ C 
ae | 
Xx 
en ee 
a 2 B+C (Wy + W 


2) 
Se esi onan eee ae eee 
ey ee ee te 


Solving the quadratic for a yields 


- ZL iL 
y, = Be (Bt do [Oy - Bi + 4)? | 
a 2c | 


But from equation (10) | 


X, = Qa + mt 
Btw). 


d 

eee x, se a = (X « = a 

" VA, x) Xs (x, x) OL BW, ) 
Let 

(X,-X -a)ad 


28 


and 


ab 
aoe (Bo - 4¢ [-cuy - Bu, + [d= B8(Wo) J] \" 


(We), ii Zc 


Grouping constants, 


ated 
, ap ae 3 
where 
ee BAe se 108 
Bee el xO 7 
r= 2295 & nee 
ee 4 
(Wo), = Ny = Hee: x 10 


The iterative process thus described results in the sequence 


mh 


| ea 
Oo a , 
a 2 _py ye 
Ky = A- (B oe 
= = . 2 
h,, A- (B CX, ) 


X 
z 


(PGE } 


wv 


A - (B= CX 
( aon 


Calculating the first few values results in a monotonic non- 
increasing sequence of points as in Figure 5. In order to 
show that such a sequence converges, we refer to the theory 
of contraction mappings as illustrated by Lyusternik and 


ids 


Yanpolskii~ 


2g 


es ss x) Xo 
Figure 5 


Lev erat cour snes, GCecralon at Lon on into ua? and 


Y= f (X). In order to solve for the equation KX =apuee 





We set up the iterative sequence of elements: 0:30 





wnere Xo iS Suwa rorurcrvary Clenedunct Bae Here we have 
Xap =t (xX ), If the sequence o.e. is convergent to some 
of 1 m ue ns 
Xx", then X* is a solution of the equation X=f (xX), and x” 
is called a fixed point of the transformation Y=f (A). 
Tne principle of contraction mappings provides a cons 
dition for the existence of a fixed point of the transfor- | 


Maud Or ey c= i ate) 


DarINITIONS: 
(1) An operator f from i into z is said to be a 
contraction mapping if there exists a constant q,(0<q<1l), 


such that for any x, x, CF, 


| f(x, ) - f(x) dl | x, -x : 


(2) A sequence ix | having the property tnat, given 
any €>0, there exists a subscript N such that |x -< | <2 fo| 
m n 


all m>N, n>N, is said to be fundamental, 
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BOLZANO-CAUCHY CRITERION: 

A necessary and sufficient condition for the sequence 
{x} to be convergent is that it be fundamental. 

The following shows tnat the principle of contraction 
mappings provides for a fixed point of the transformation 
mee t (Xs 

Lev f be a continuous operator from E into ae such 
that f is a contraction mapping and in an 1.€. twice con- 
tinuously differentiable in all components. MThnen there 
exists a solution X” to the equation <= T(x), and this 
Solleron is unique. The ltverative sequence formed by 


successive approximations is convergent to X* whatever the 


initial approximation Ko- 


BROOr : 


Form the iterative sequence of elements 


- = W 

4 = F(X) 

ky = £(X,) 

kK 2 £(K, 4) (11) 
Koei = Fh) (12) 


Subtract equation (11) from (12) and take absolute values 


Kar 7 Fabel tO) - 2S) 


of 


But f is a contraction mapping by hypothesis. 


i 


| An 4d ‘ x | 


pe(4)- £1) (24]% ~ ee 


rs 
[A 


les eee ee (13)| 


eee ~ x | 





2 
KOR SG Kl | ear aie 
ey ee a Z ow aL O 
Nee ~x lequix -x | (14)) 
n+ 1 he L O 


If the sequence, equation (14), can be shown to be funda- 
| 
mental, tnen by the Bolzano-~Cauchy criterion it is convergel 


Choose €>0, and N large. To show that | as i < ¢ for all 
n | 


m > N, n>N, we add and subtract all values of X,, m<k<n,| 


Xx -~ X + XK ~ X of X ae ee "FE Ga = ] 


< |X =X 


; | + | x 
Mm i= m=] 


am 


a a * © 9 “f- AK oa x 
m-2 ne n 





From equation (14), 





| m= 1 fe { 
|S ga | 
M=2 4. ] 

| mi ~4,2.|24 | 4° 4, | 


x = A 


am | fox 
Nef] 7 Soe | ioe 


De 


[Xa gh ge Xp gt K yg Kol Sara bes gy git )s [xq “x 


m<-1l m-l meses 7 


Ing-n=1 


Yh 
< a (l+qra@+> ora eee 





Seance q <= 1, the partial geometric series is less than 
ee 


nal 1 
-~ X of > = 000 -= ° cs i ee a 
m ta~L ea en el ‘| 4 nL 


ee 





Ve can choose n large enough so that the sequence is funda- 
mental and therefore convergent. Since q<1, we can make 

the right hand side arbitrarily small. In fact, | 4 - 4,.4h¢9 
aS N-Pro, 


From equation (11) 


= P(x ) 


1 
| 4.7 2(K)) | =|f(x ,)-f(4 ) <a x - x 
[x,- £(<,)| <a[k,-<4 [> 0 

Since f is in 0°, as X—r x", f(x,)f(4*). 


iP anes: 


| eae) 





and 


x* = £(X*) 
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if ere <“q<1, then f is a contraction mapping. From 


the mean Value theorem of differential caiculus 





f(b)-f(a) = (b=a) « £*(X) for some X, a <= tae 
cr (gy = £62) = fla) 
bea 


if | frt(x) | <q <1, then|f(b)-f(a) |< |o-al 
and | f(b)-f(a) |< q|o-a] | 


From equation 10,1 the function under study is; 


toll 


(f)= A - (B- CX) 


ieee? | = eer 
2(B=-CxX)* 





To show that the function is a contraction mapping it is 
necessary to show that [f'| < q < 1 for all values of X 


that will be encountered, It is easily seen that when 


A = - the derivative is not defined, and tnat when X > B 
C 


it is imaginary. 


The values of the constants in the equation are; 


A = 5.65 x 107 
B = 1.81 = 107 
C = 2.35 x 10% 
Kad, = 1.6 x 10% 
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HesGecilre ver cnow that <a * = (ae ae & 10° for all X that 


are encountered in the iterative process, The first 


approximation to X is; 


Vas) ae 4 
Ko = iE = 1.6 x 10 


= 5,65. x Or (ie Clee ho? oe 235 102 eee 10 )e 


x, = 5,65 x 10' - (18,3 x 108 )z 


4 


= Ly 
x = 1,39 x 10° < Xo 


SRGGe iO Teele o Ler 10! ae, 35) 10S 5a 


i 
2 


Ps 
a) 
| 


i) 


a uy 
i oe Oe = x < Xo 


rs 
| 


It can easily be shown that each value of X will be smaller 
than the one that precedes it, since at each iteration a 
larger value will be subtracted from the constant A. Thus 
the sequence Ay is monotonic non=-increasing, and its upper 
bound is W,. Since ui, <-e , all values of X that will be 
encountered are also strictly less than : 2 Solving for 
the first derivative when X = q, yields: 

2.35x10° 


io 
2(4,243 x10") 


eg! |= 
Thus | f° | < 1 for all values of X in the iterative sequence. 
Mis is so because { Xk is monotonic non-increasing, and 
any value of X < Ny increases the denominator of the deriv- 
ative, thereby decreasing its value. Thus f is a contrac-~ 
tion mapping, and the sequence ix} converges. Lyusternik 
and vepeleaae show that if f is a contraction mapping, 


then the sequence is convergent to X* as fast as a geometric 


oe. 


progression with ratio q. Thus the computer iteration 
procedure is not lengthy since in this case q = .0022 

The imowledge that we can solve for X,, and thus Xs 
at the start of the horizontal leg of the flight assures us 
that we can correctly find the average weight on the leg 
using the procedures described earlier. We can solve 
Wy - Fs Woe Knowing Wo we solve for the fuel used in the 
descent using F= f (Wo9h), and N3 =Wo- F. The fuel 
burned during this profile CW. - a3) is stored, along with 
the altitude and proper airspeeds for the profile. 

Once fuel figures for all profiles are calculeved tae 


searching procedure selects the global minimum and prints 


it as the optimal fuel along witn corresponding altitude 





and airspeed figures. The possibility of a tie is Virvuem 
eliminated by reading fuel figures to four declines 2¢--y 
should several local minima arise, the search is desmamed 
so that the global minimum is always selected, 

A listing of the Fortran IV program and a sample out- 
put is given in Appendix B. All variable names are written 
such that they are easily recognizable, TAS is true airspe 
DIST is distance, etc. As inputs, the program only require’ 
the course and distance to the field, the type aircraft 
(A4 = 1, A7 = 2, etc, This allows for the difference in 
aircraft), the weight of fuel aboard the aircraft at time o 
BINGO, and the wind information. The wind can be up-dated 


at intervals as set by meterological readings. 


CHAP hs ti 
RESULTS AND CONCLUSIONS 


For a fixed distance, X, nautical miles, (100 < XK, < 250), 
the winds aloft and the aircraft weight combine to determine 
the optimal flight profile. Under a zero wind condition 
and normal BINGO fuel weight, the optimal flight profile is 
to climb to the service ceiling and commence an idle descent 


pe Or CO Clmmbp Until xX. is reached, and then to start 


d. 
the idle descent such that xX is zero, For fuel weights 
higher than the normal BINGO weights, the optimal altitude 

ms Lower than the service ceiling because of the excess 

fuel used in climbing to high altitudes when the aircraft 

is heavy. 

Since zero Wind conditions are rarely if ever en- 
countered, the wind velocity aloft is a most important factor 
iieactcermining the optimal flight profile. In order to 
test the effect of winds on a standard BINGO problem, a 
typical wind pattern for winter months along the coastal 
region of Southern California was obtained from the weather 
facility at the Naval Auxiliary Landing Field in Monterey. 


The following winds are typical of a northern hemis- 


Ppaere CYCiOnNIG Low pressure area. 


EN 


OGIO R ULES WIND 

G2] 5000 tt 270°/ 10 eee 

6 =10000Ntt 24U0°/7 15 kts 
ay =e OCe rs 210°/ 20 kts 
16 =Z20000Mrt 1809/7 30 ite 
21 ~25000 ft 150°/ 40 kts 
26 = 30000 1 t 1209/7 60 kes 
31 -35000 ft 095°/ 80 kts 
36 -40000 ft 080°/ 100 kts 


The BINGO situation considered is such that a pilot Leaves | 
a carrier off the coast of Southern California and flies 
090° / 160 miles to NAS North Island near San DIEGO wate | 
fuel weight at start of BINGO is 22004. 

For this example a flight profile with the service 
ceiling as the optimum altitude would require 1499.4#%. A 
guess by a "seat of the pants" acquaintance of mine would 
use 25000 ft as the optimum altitude and would require 
1435. 3#- 
flignt profile as follows, 


The computer flight plan predicts the optimal 


Optimum altitude is 
Climb speed is 
Cruise speed is 
Start descent wnen 
Descent speed is 
Fuel yrequired is 


The savings in jet 


the guess, and 364% for the optimal over the climb to serv: 
ceiling is substantial when it is considered that this 
Savings would allow two more landing attempts either at thi 


carrier or at the field ashore, 


314.4 
250.0 kts 


35:5" 


Lingo LU 
kts 
Set aCo OUT 


kts 
lbs 


fuel of 300% for the optimal over 








Multiply this savings by 
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thousands of flights per month if the system were incorpor- 
ated in the fleet, and it is easily seen that the jet fuel 
saved would soon become Significant. The two extra landing 
attempts would probably resuit in fewer aircraft ramaining 
on the beach over night, 

In order to demonstrate the usefuiness of extending 
the present BINGO computer program to long distance flights, 
a flight of 600 miles was flown with a fuel weight of 6000 
pounds at the start. The entire distance is flown with a 
course of 115° magnetic. This is not realistic since even 
West - Hast cross-country flights sometimes require course 
changes of more than 20°, Also, the present program does 
not change altitudes to nigher ones when the aircraft weight 
decreases to a value that suggests a climb in order to 
remain optimal, Nevertheless, some fuel values for 
different profiles indicate the savings in fuel that would 


Pemrecaiezed it the Optimal profile is flown, For example; 


Fuel (lbs. ) Altitude (ft. ) 
5596.6 5,000 
4646.1 10,000 
4333.9 15,000 
-—>~ 3842, 3 20,000 
4517.1 25,000 
4346.8 30,000 
3976.6 35,000 
oa 4O,000 Service 
Cer line 


a 


Im this case the service ceiling is Che optimal 
altitude, But if this altitude can not be reached because 
of poor engine performance in a particular airplane, then 
20,000 ft. would be the best altitude, 

An experienced pilot may often fly a near-optimal 
flight profile merely by an educated guess or by carefully 
Managing the fuel flow, But he can never do better than 
the computer prediction, and will probably do worse most 
of the time since his wind information is sketchy if he 
has any at all. 

The professional doubter may complain that even mete 
Ology doesn't know the accurate wind information at altitu 
so why even bother with a computer solution? We can only 
reply that some good estimate is better than none at all, 
and that modern meteorological equipment can measure and 
predict the winds quite accurately. 

The incorporation of a computer system to predict 
optimal flignt profiles may at first meet with some inerti. 
from fleet pilots, especially the more experienced ones. 
The “seat of the pants" pilot may reject the computer 
decision as hocus-pocus or just plain incorrect. But even 
tne most experienced aviator may encounter vertigo some 
dark, rainy night and he may discover that flying the 
airplane is about the only job he can handle, Figuring ou 
an optimal flight profile to the beach may take a back 
seat to survival, It is in this situation that informatio 


from departure control would be most welcome, 
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The cost involved in establishing and operating a 
System such as described above may appear to be too large 
when compared to the small savings in relatively inexpensive 
jet fuel. But the possible loss of a Phantom II and/or 
a pilot may shift the balance in the system's favor.) 

One possible reaction against the use of computer 
flight planning for cross-country flights is that often the 
air traffic control center will not allow a climb to the 
desired altitude because they must wait until they can safely 
fit the aircraft into the West-East (or vice-versa) traffic 
flow, which may be considerable in this jet age. This will 
result in remaining at a non-optimal altitude for many 
minutes, But this drawback exists for non-computer flights 
as well. And it is still better to climb from a non-optimal 
altitude to an optimal one wnen cleared to climb, than to 
go from one non-optimal altitude to another, 

Another great advantage in long range computer flight 
planning occurs wnen several alternate routes are available. 
Then the program can be written to choose the optimal 
altitude as well as the best route of flight. The many 
possible routes and altitudes would require too much compu- 
Maelon ror a pilot during his flight briefing, but a high 
speed computer can solve such a problem in minutes, Thus 
the pilot can spend nis pre-flight time briefing the mission 


and tactics, and let the machine do the arithmetic. 
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Aircraft Performance Charts for the Navy ASC Aircraft 
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MODEL: A-4A, A-4B, A-4C DATA AS OF: 1 DECEMBER 1965 
ENGINE: J65-W-16A DATA BASIS: FLIGHT TEST (NAVY) 
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Figure 6. Climb Fuel 
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PRESSUNE ALTITURE — 1000 FEET 


MODEL: A-4A, A-4B, A-4C 
ENGINE: J65-W-16A 


DATA AS OF: 1 DECEMBER 1965 
DATA BASIS: FLIGHT TEST (NAVY) 
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MODEL: A-4A, A-4B, A-4C 


ENGINE: J65-W-16A 





INITIAL GROSS WEIGHT — 1000 POUNDS 
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Figure 8. 
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DATA AS OF: 1 DECEMBER 1965 
DATA BASIS: FLIGHT TEST (NAVY! 


Climb Time 








AVERAGE GROSS WEIGHT — 1000 POUNDS 
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MODEL: A-4A, A-4B, A-4C 


ENGINE: J65-W-16A 
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DATA AS OF: 1 DECEMBER 1965 
DATA BASIS: FLIGHT TEST (NAVY) 
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CRUISE —— FUEL 


MAXIMUM RANGE 


DATA AS OF: 1 DECEMBER 1%65 


DATA BASIS: FLIGHT TEST (NAVY) 


MODEL: A-4A, A-4B, A-4C 
ENGINE: J65-W-16A 
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Figure 10. 
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GROSS WEIGHT — 1000 POUNDS 


DISTANCE — NAUTICAL MILES 
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MODEL: A-4A, A-4B, A-4C 
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DESCENT SPEED SCHEDULE — KCAS 





GROSS WEIGHT — 1000 POUNDS 


0 


180 200 
100 165 180 
200 155 170 





Figure 11, Descent Distance 
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DATA AS OF: 1 DECEMBER 1965 
ENGINE: J65-W-16A DATA BASIS: FLIGHT TEST (NAVY) 
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MODEL: A-4A, A-4B, A-4C 
ENGINE: J65-W-16A 


Figure 12. 


DATA AS OF: 1 DECEMBER 1965 
DATA BASIS: FLIGHT TEST (NAVY) 
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DATA AS OF: 1 DECEMBER 1965 
DATA BASIS: FLIGHT TEST (NAVY) 
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APPENDIX B 


Program BINGO and Sample Output 
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